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Abstract 



> 

. The finite part of the self-force on a static scalar test-charge outside a 

Schwarzschild black hole is zero. By direct construction of Hadamard's ele- 



I mentary solution, we obtain a closed- form expression for the minimally cou- 

I pled scalar field produced by a test-charge held fixed in Schwarzschild space- 

■ time. Using the closed-form expression, we compute the necessary external 
O^' force required to hold the charge stationary. Although the energy associated 
$H . with the scalar field contributes to the renormalized mass of the particle (and 

thereby its weight), we find there is no additional self-force acting on the 
charge. This result is unlike the analogous electrostatic result, where, after a 
^ ■ similar mass renormalization, there remains a finite repulsive self-force acting 

■ on a static electric test-charge outside a Schwarzschild black hole. We confirm 
our force calculation using Carter's mass-variation theorem for black holes. 
The primary motivation for this calculation is to develop techniques and for- 
malism for computing all forces — dissipative and non-dissipative — acting 
on charges and masses moving in a black-hole spacetime. In the Appendix we 
recap the derivation of the closed-form electrostatic potential. We also show 
how the closed-form expressions for the fields are related to the infinite series 
solutions. 
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I. INTRODUCTION 



In order to gain deeper quantitative understanding of highly relativistic binary star 
systems spirahng toward coalescence, a number of authors have used perturbation theory 
to study the motion of test-particles orbiting black holes. (See e.g. [0-^.) The basic idea 
is to solve the linearized field equations with out-going radiation boundary conditions on a 
Schwarzschild or Kerr black-hole background. The source of the perturbing field is generally 
chosen to be a test-particle (endowed with a small scalar charge, electric charge, or mass 
charge) moving on a bound geodesic orbit, e.g., a circle. Once the perturbing field (scalar, 
electromagnetic, or gravitational) is calculated, the energy radiated per orbit {i.e., the time- 
averaged energy flux) can be computed by a performing a surface integral over a distant 
sphere surrounding the system. The rate energy is carried away by the radiation is then 
equated to a loss of orbital energy of the particle; thus one can compute the rate of orbital 
decay.Q In effect, this energy-balance argument gives the time- averaged radiation-reaction 
force associated with imposing out-going radiation boundary conditions 0. This method 
has been used successfully to compute the inspiral rate of coalescing binaries to very high 
relativistic order, and it has been used to check analogous post-Newtonian calculation of 
the inspiral 0. These perturbation calculations also add insight into effects such as wave 
"tails". (See |]T| and |^ for post-Newtonian and perturbation-theory discussions of tails.) 
However, in spite of the success of the perturbation approach, the method, as it has been 
applied, has a drawback: it has only been used to compute time- averaged, dissipative forces 
on the particle. 

The time-averaging is problematic for two reasons: (1) During the late stages of inspiral 
the orbit will be decaying swiftly, and there will not be many orbits left before the final splat 
to average over. (2) In the curved spacetime near a black hole, the self-forces experienced 
by the particle are not instantaneous forces. Rather these forces arise because the fields 
produced by the particle at one instant travel away from the particle, encounter the curvature 
of the spacetime, and then scatter back and interact with the particle at a later time. (See 
DeWitt and DeWitt p for an illuminating discussion of this point.) Thus the fields produced 
when the particle moves through, say, periastron, and suffers the maximum coordinate 
acceleration, will not come back to interact with the particle until later in the orbit when 
the coordinate acceleration will be more gentle. In other words, the state of motion that 
produced the fields, will not be the state of motion affected by the fields. Although it is 
unlikely that taking into account the non-local nature of the self-forces will qualitatively 
change any of the results found by the "averaging" technique, the issue has never been 
rigorously addressed with perturbation theory. 

In addition to the dissipative forces, the perturbing fields will produce other non- 
dissipative — conservative — self-forces on the particle, and these forces have not been 
studied by the conventional perturbation-theory approach [As DeWitt and DeWitt 

have demonstrated, these forces can be computed using a different perturbative approach.] 
For example, suppose we have an electric charge in a circular orbit around a Schwarzschild 



Finding a gravitational- wave signal in noisy detector data requires an accurate prediction of the 
orbital decay rate. See Thorne p. 
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black hole. Even if we neglect the dissipative effects of the radiation reaction, which cause 
a secular decay of the orbit, the charged particle still does not travel on an exact geodesic 
of the spacetime. The particle will feel an additional conservative force — proportional 
to square of the charge — that pushes it slightly off the geodesic. As this force is not 
radiative in nature, the portion of the field that produces it is explicitly discarded in the 
conventional perturbation approach because it falls off faster than 0[l/r], and therefore any 
effect it might have on the motion is missed by a technique that uses a surface integral on a 
distant sphere to determine the back-reaction effect. (2) Another perturbation calculation 



by Gal'tsov ITQ] has used the "half-retarded-minus-half-advanced" Green's function to find 
the radiation reaction forces. But these calculations also miss the conservative {i.e., time 
symmetric) parts of the force, because they simply cancel when the two parts of the Green's 
function are subtracted. In other words, the Gal'tsov result differs from the DeWitt-DeWitt 
result in that it misses the conservative part of the force. 

We can try (but as we will see, fail) to understand the origin of the conservative forces by 
considering a static electric test-charge |Tl|] held fixed by some non-conducting mechanical 
struts outside a Schwarzschild black hole. The charge experiences a self-force (which is 
clearly not a radiation-reaction force in any conventional sense) proportional to the square 
of the electric charge. A physical origin of such a force might be naively described as follows: 
the surface of the uncharged black hole will act as a conductor, and therefore the presence 



of the external charge will induce a dipolar charge distribution on the horizon [0. The 
magnitude of the induced dipole moment will depend on the size of the dipole (in this case, 
the mass M of the Schwarzschild black hole) and the magnitude e of the external charge. 
This leads one to believe that the charge will experience an attractive dipolar force scaling 
as b~^, where bg is the Schwarzschild radial coordinate of the charge. This force will be 
in addition to the natural attractive force (the weight of the particle) which, according to 
Newton, should scale as 67^. Thus our naive suspicion is that, in addition to the weight, the 
strut holding the charge fixed will have to counteract an attractive self-force of the form 

^(naive) ~ "7^ • (1.1) 

Smith and Will (SW) [|l3l have calculated the force required to hold an electric charge 
fixed outside a Schwarzschild black hole. They find the total force exerted by the mechanical 
strut holding the charge fixed must be 



^(strut) = ^^(^l-^j (1.2) 



Here /iren is the renormalized mass of the test-charge. We use units in which G = c = 1. As 
expected, the first term scales as and shows that the strut must support the weight of 

the particle; such a term would be present whether or not the particle is electrically charged. 
The second term shows — contrary to the "physical intuition" presented above — that the 
black hole tries to repel the charged particle and the repulsion scales as 0[67^]. SW show 



Although the second term in Eq. (1.2) appears only to be accurate to leading order in M, 



for almost any choice of variables in Eq. ( |1.2| ) the first term dominates, and the second term 
only lightens the load a small amount. 

The fact that our naive physical intuition led to both the wrong direction and the wrong 
scaling of the force suggests such calculations should be carried out in detail before conclu- 
sions are drawn. In this paper we carry out the analogous calculation for a scalar charge 
in the presence of a Schwarzschild black hole, i.e. we compute the force required to hold a 
sca/ar charge fixed. The primary result of this paper is to show that, although the scalar field 
does contribute to the renormalized mass, there is no force analogous to the second term in 



Eq. (1.2) for a scalar charged particle. The primary motivation for this calculation is much 



broader: to begin to develop techniques and formalism for tackling a sequence of problems, 
namely, calculating all the forces, linear in the fields (dissipative and non-dissipative), that 
act on test-charges and masses moving in the proximity of a black hole. 

As the implicit underlying motivation for this work is calculating the forces on moving 
charges, we can ask: does our no-self-force-on-a-statzc-scalar-charge result carry over to 
the case of a moving charge? Not directly. However, we can gain some intuition about 
extending our static results to dynamic results for scalar charges by examining the static 
and dynamic cases for electric charges. Although the SW force calculation is only valid 
in the static limit, DeWitt and DeWitt P] show that when the electric charge is in slow 
motion around the black hole, in addition to the conventional radiation reaction force, there 
is a conservative repulsive force acting on the charge. This force is independent of the 
velocity, and clearly corresponds to the force found in the static calculation, i.e., the second 
term in Eq. ( |1.2|) .P| Therefore, it seems reasonable to expect a correspondence between the 
static case and the slow-motion case for scalar charges. For scalar charges, this means the 
absence of a self-force on a static charge should imply there is no conservative (repulsive or 
attractive) force on a slowly orbiting charge. This, however, does not rule out the presence 
of a higher-order, conservative force that scales as, say, the square of the orbital velocity. 
(This issue is under vigorous investigation [l^.) As a consequence, although both electric 
and scalar charges will spiral inward due to radiation reaction forces, scalar charges will not 
suffer the same persistent, velocity- independent, conservative force trying to hold them off 
the geodesic. 

As a matter of principle, many of the delicate issues of radiation reaction forces in 
curved spacetime we will be dealing with in this paper have already been addressed in quite 
a general context by DeWitt and Brehme [0 (using a world-tube calculation), by Quinn and 
Wald (using an axiomatic approach), and Mino, Sasaki and Tanaka [jrU| (using matched 
asymptotic expansions). However, in practice, the results presented in these papers have 



never been applied to give even simple results like Eq. (1.2), nor did they work out the 



the result is exact for static charges; the appearance is only a beautiful artifact of Schwarzschild 
coordinates. Switching to harmonic or isotropic coordinates p^JT5|| makes the second term appear 
as an infinite series in M. See Eq. (|]|b). 

^ DeWitt and DeWitt only compute the force to leading order in the mass of the hole and their 
result appears to agree exactly with the second term in Eq. (p..2|). However, the agreement is really 
only at leading order in M. See the previous footnote. 



A 



general expressions for forces on scalar charges. 

We begin in Sec. |I| by calculating the scalar field produced by a static scalar charge in 
Schwarzschild spacetime. The standard method for solving the field equations on a black-hole 
background is to use separation of variables and decompose the solution into Fourier, radial 
and angular modes (spherical harmonics). (See Appendix ^ for an example.) The result is 
an infinite-series solution for the field. This method works well for computing the far-zone 
properties of the field {e.g., the energy flux) where the field is weak and can be accurately 
described by the first few terms in the series. However, when computing a force, such as 
Eq. (|1.2| ) (or its scalar analogue, as we are doing here), the behavior of the field on a distant 
sphere will not suffice; we need to know the detailed behavior of the field up close to the 
particle. For this it will not be sufficient to describe the field by a few multipoles. Describing 
the singular nature of the field near the charge requires an infinite number of multipoles 
20| , |2l[| ; therefore it is not surprising that SW did not use the standard multipolar expansion 



of the electrostatic field {e.g., p2| , pB| ) in their force calculation. Rather, when embarking on 
their electrostatic force calculation SW comment: "... an exact calculation is made possible 
by the fortuitous existence of a previously discovered analytic solution to the [static] curved- 
space Maxwell's equations" . The exact solution to which they refer is an old result due to 
Copson and modified by Linet ||2^ which gives a closed-form expression for the electrostatic 
potential of a fixed charge residing in Schwarzschild spacetime. [See Eq. ([B5|).] The Copson- 
Linet formula has the advantage that the singular nature of the field near the electric charge 
is manifest, thus allowing simple calculations in the close proximity of the particle. We 
begin by constructing the scalar- field analogue of the Copson-Linet result: a closed-form 
expression for the field of a fixed scalar charge in Schwarzschild spacetime. Our derivation is 



similar to Copson's, and is based on constructing the Hadamard "elementary solution" ||25 
of the scalar-static field equation. As in the Copson-Linet electrostatic solution, our solution 
will clearly show the divergent behavior of the field near the particle. In Appendix 0, we 
outline Copson's derivation of the closed-form expression for the electrostatic field with the 
Hadamard formalism. In Appendix M, we generate some interesting summation formulas by 



equating the closed-form solutions with the infinite series solutions [22,23 



As a historical note, we mention that Copson obtained his solution for a static 
electric charge in a Schwarzschild spacetime in 1928. This is approximately 40 years before 
the term "black hole" was coined [^ and researchers began formulating the no-/iazr theorem. 
However, Copson was probably the first to note one of the important features of the no-hair 
theorem: "... the potential of an electron on the boundary sphere R = a [on the horizon] is 
independent of its position on the sphere, a rather curious result." 

Our force calculation in Sec. HI closely follows the SW calculation of the electrostatic 
force. We use the exact expression for the scalar-static field found in Sec. ^ to compute the 
stress-energy tensor and the force density in a local, freely-falling frame near the fixed scalar 
charge. By integrating the force density over a small spherical volume of radius e centered 
on the charge (spherical in the freely-falling frame) we compute the total force the strut 
must supply to hold the charge fixed. In the limit where the radius of the spherical region 
of integration shrinks to zero, there remains a formally divergent piece of the force scaling 
as (charge)^/e. This factor multiplies an acceleration that is identical to the acceleration in 
the first term in Eq. ( |1.2|) ; therefore the infinite term is simply absorbed into the definition 
of the mass of the particle, i.e., in Eq. (|3.30|) we define 



1 g2 

yWren = /^bare + - lim — , (1.3) 

where q is the scalar charge of the particle, and /ibare is the bare mass of the particle. In this 
way, the scalar charge does contribute to the renormalized mass of the particle in exactly 
the same way as the electric charge contributes to the renormalized mass in Eq. (|1.2|) . The 
form of classical mass renormalization depicted in Eq. (|1.3|) is seen in all calculations of this 



type {e.g., [^,^,0). The main conclusion of this paper is that, after the formally infinite 
piece is absorbed into the renormalized mass, there remains no scalar counterpart to the 
second term in Eq. ( p.. 21 ). Also in Sec. |IT|, we verify this answer using conservation of energy. 

Neither the charge's contribution to the renormalized-mass term or to the repulsive term 
in Eq. (|1.2|) should be confused with the contribution to the gravitational force on the 
particle due the stress energy of the electric field perturbing the metric of the spacetime. 
Such a force would scale as /xe^. In this calculation we are explicitly ignoring corrections to 
the metric. See the discussion following Eq. 



In Sec. 1^ we discuss a number of alternative methods for solving similar problems, 
as well as similarities and differences in the scalar-static and electrostatic results. We 
also collect what is known about the forces on static charges (mass, electric and scalar) 
in Schwarzschild spacetime. This gives a clear indication of what future research is needed. 



II. SOLUTION OF THE SCALAR FIELD FOR A FIXED CHARGE IN 
SCHWARZSCHILD SPACETIME 

In order to find the self-force acting on the scalar charge, we first solve the massless 
scalar field equation 

DV ^ (l/^)(^g-^V^)^^ = Anp (2.1) 
in Schwarzschild spacetime.0 Here commas denote partial differentiation, Greek indices run 



to 3, and Latin indices will run 1 to 3. The source p for our field will be a point-like pS 
scalar charge which can be described by 



(2.2) 



p(t,x) = q / (l/V^) 5\x--h'^{T))dT 

J —oo 

q 53(x-b(t)) 



where 6"(r) is the spacetime trajectory of the charged body, r is the proper time measured 
along the path, andw* = dlP/dr = dt/dr. [Notice, / p(i(proper volume) = / pu^^—gd^x = q 



^Since the Ricci scalar curvature R vanishes in the Schwarzschild spacetime, it would seem that 
including coupling to the curvature (e.g. a conformally invariant term ^RV) in Eq. ( p.lj ) would 
have no effect on our results. However, if we include coupling to the curvature, the stress-energy 



that enters the force calculation [Eq. (3.9a)] would have to be modified also. Therefore, we make 



no claim that our results hold for non minimally coupled fields. 



in a frame comoving with the charge.] We will later restrict our attention to a static field 
and a fixed source charge at b = 6z, but, for the present, we will leave the time-dependence 
in the equations. 

We use isotropic coordinates [jl4| to describe the Schwarzschild geometry. The line ele- 
ment is 

In these coordinates Eq. can be written 

C'^Vij{uJ, x) + C^Vjiuj, x) + CV{uj, x) = 47r(l + M/2r)^p(cu, x) , (2.4) 

where 

C'^ = diag(l, 1, 1) , (2.5a) 



= hir)— = — ln(l - iM/2rf) — , (2.5b) 



, (l + M/2r)6 

^ = " (1 - M/2rr ' ('-'^^ 



and we have used the Fourier transform 

V{uj,^)e-''^'dLJ (2.6) 

-oo 

to eliminate the time variable. We can also write Eq. ( |2.4| ) in the more compact form 

V^V + h{r)Vr + C(r, uj)V = Ait{1 + M/2r)^p{uj, x) , (2.4') 

where is the fiat-space Laplacian. Our primary attention will be focused on the static 
— zero frequency — case (tu = C = 0), but the formalism we are developing is valid for the 
"Helmholtz"-type equation in Eqs. ( p.4|) and ( |2.4f ). Following Hadamard (p. 92-107), the 
elementary solution [i.e. the Green's function [^]) for Eq. (p.4|) takes the form 



Uelen. = T'^'^U^ + U^T + U^T^ +...), (2.7) 

where F is the square of the geodesic distance (in the sense of the "metric" C^^ = 5^^) from 
the source point x' to the field point x, i.e., 

T = {x-x'f + {y-y'f + {z-z'f. (2.8) 

Since we are working in three dimensions (an odd number of dimensions), there is no natural- 
log term (InF) in the elementary solution. The f/ri's are non-singular functions everywhere 
outside the horizon. Recall, in isotropic coordinates |T^ the horizon is at r = M/2. The 
simple form of F in Eq. ( p.8| ) is a consequence of the isotropic coordinates we are using; 
therefore the technique we are developing cannot be easily extended to geometries that do 
not have spatially isotropic coordinates, e.g. Kerr geometry. 

The formula for the leading order behavior of the series is given by Hadamard (p. 94) 
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?7o(x, x') 



= exp 



C'^T a + c-^r , - 6 



dX' 



(2.9) 



where A is the arc length measured along the geodesic connecting the source and field points 
Letting 6 denote the angle between the two spatial vectors x and x' and using Eqs .( p.5| 
and (2^) we have 



t/o(x,x ) = exp I - [h{r){r - r'cos6')] — \ 
= exp|-£Mr")^} , 



1 - (M/2r')- 
\\ 1 - (M/2r)2 



(2.10a) 
(2.10b) 

(2.10c) 



where r = |x|, r' = |x'|, and r" is the dummy integration variable over r. We have also used 
the geometric relationship (r — r' cos 6)dX = Xdr. 

Because we are looking for an axially symmetric solution, we may assume the f/„'s are 
functions of only the radial variables [e.g., notice Uq = Uo{r, r')]. We now substitute Eq. ( [^.7[ ) 
into Eq. ( p.4| ), use relations such as ^^-^T^jFj = 4r, and collect powers of T. The result is 



i2Uo,r + hUo) 



^2 lyt^'^ 

2r 



n=0 



(2n+ l)(2t/, 



n+l,r 



n+l) 



^2 _ 

2r 



1 2ri — 1 
+2(2^2 + n - l)t/„+i + -{rUn),rr + /it/„,r + (2t/„,. + /it/„) + C(r, oj)[/„ 



■pn-l/2 _ 

(2.11) 



Since Ueiem is a solution to the homogeneous equation everywhere (except at the source 
point) the coefficients of each power of F must vanish independently. The first term gives 
an equation for Uq 



2Uo r + hUo = 



(2.12) 



Notice that our Eq. ( p^) has already given us a particular solution (Eq. (|2.10c| )) to this 
equation. Setting the coefficient of r"~^/^ to zero and using an integrating factor, we obtain 
a recursion relation for the f/„'s 



Un+ijr, r') 
Uo 



-1 



2 _ ^/2\n+l 



2n + 1 



{T'"Un) r"r" 2r"UQ r"Un r"r" 



+ i2n~l)iUjUo),r"+r"Cir",uj){UjUo) 



dr" . (2.13) 



In the integrand the UnS are functions of the dummy integration variable r" and the source 
point r', i.e. Un = Un{r",r'). This recursion relation allows us to construct a series solution 
for the Green's function for the Schwarzschild-Helmholtz equation Eq. ( p.4|) or ( p.4f ).p| 



We can gain confidence in the recursion relation Eq. ( 2.13| ) by applying it to the true Helmholtz 



S 



We now explicitly assume that the charge and the field are static, i.e. we assume that 
to" = C(r, cu) = in Eq. ( p.l3|) . Beginning with Uq from Eq. (|2.10cj ) we can construct 





1 


(2M)2 




Uo 


2 (4r'=^ - 


M2)(4r2 - 


M2) ~ 


U2 


3 


(2M)4 




Uo 


8 (4r'2 - 


M2)2(4r2 - 


-M2)2 


U3 


5 


(2M)6 




Uo 


16 (4r'2 - 


-M2)3(4r2 


-M2)3 



'2^ 
3 

o7 



16 
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(2.14a) 
(2.14b) 
(2.14c) 



Substituting these into Eq. (p.7|), the pattern is immediately clear and the summation is 
elementary 



U^( 
Uo 



v/rv/TTT> 

1 r(4r'^ 



M2 



r' J(4r'2 - M2)(4r2 - M2) + 4M2r 



(2.15a) 
(2.15b) 

(2.15c) 



Equation (|2.15c| ) is a closed-form expression for the Green's function for Eq. (|2.4| ) (with 
C=0). In the limit M = it reduces to |x — x'|^^, the Green's function for Poisson's 
equation in fiat space. We use Eq. (|2.15c| ) in precisely the same manner we would use 
the Green's function for Poisson's equation: we integrate it against the source to obtain a 
particular solution to the inhomogeneous equation. In our case the source is the static scalar 
test-charge held fixed at the spatial position b, and we have 



l^art(x, b) 



-(1- 



2b -M 1 



2b J 



Uelem{^,^)d X 

4rb 



2b + M /r(x, b) 7(462 - M2)(4r2 - M^) + 4M2r(x, b) 



(2.16a) 
(2.16b) 



equation in flat spacetime (Eq. (2.4) with M = 0, but C = lu'^ ^ 0). In this case, Eq. (2.9) 
gives Uo = 1, and the recursion relation Eq. ( 2.13 ) gives Ui = —uj'^/2, U2 = a;^/24, C/„ = 
(— l)"u;2"/(2n!). Summing the series gives the Green's function for the Helmholtz equation 



U, 



elem 



cos(a;|x — x'| 



Now compute the inverse Fourier transform and we have 



G{x, x ) 



1 I 6[t' - (t- |x-x^|)] ^ 5[t' - (t + |x - x^l 



X' 



X — X' 



This is the half-advanced plus half-retarded Green's function. See |21|, Eq. (6.61). 



Q 




2rhhh cose + hi -M^ sm^ 9 



(2.16c) 



In the last two steps we have exphcitly included the factor l/u^{b) = y—goo{b) = {2b — 

M ) / (26 + M) = ^ iph — M) I {bh + M) ; in the last step we have converted to harmonic 
coordinates [|I4|. The leading minus sign in Eq. (|2.16|) is a consequence of the source term 



in Eq. (|2.1| ); our source is '+'47rp and not the familiar '— '47rp of electrostatics. (Like scalar 
charges attract.) As we are now dealing strictly with a static solution, we have also dropped 
the twiddle denoting the Fourier transform in Eq. (|2.4| ). 

Equation ( |2.16| ) is a particular solution to the scalar-static field equation, but is it the 



desired solution to the equation? In other words, does it satisfy all the boundary conditions? 
First, the field and its derivatives are well behaved outside — and on — the horizon; thus 
our solution has no unphysical regions of infinite energy (save, of course, at the location of 
the charge). For r ^ b > M/2 we see from Eq. (p.l6[ b) 



W--oo)--^r^^). (2.17) 



.26 + M. 

As the charge is lowered toward the horizon, the factor (26 — M) / {2b + M) extinguishes 
the field measured by a distant observer. Notice if the charge is lowered to the horizon 
{b = M/2) the field completely disappears. (See, e.g., for discussion.) Thus, it is the 
extinction factor — which had its origin in the factor of l/u^{t) in Eq. ( |2.2D and has survived 
throughout the calculation — that enforces no scalar hair on the black hole^ 

Equation (|2.17|) is the appropriate asymptotic form of the scalar field; therefore the par- 
ticular solution Eq. ( |2.16| ) is the desired solution which satisfies all the boundary conditions 



32,33 



V{x,h) = Vpart{^,h) . (2.18) 

Before embarking on the force calculation, we note a remarkable feature of the closed- 
form solution we have found. Despite the fact that every term in the Hadamard series is 
divergent at the horizon [Eq. (|2.14|) ], the closed- form expression for the elementary solution 
is well behaved on the horizon. The easiest way to see this is with the harmonic-coordinate 
expression Eq. (|2.16| c) evaluated on the horizon {r^ = M) 



It is also interesting to note that the horizon is not a surface of constant "potential". In 
Sec. I our (flawed) intuition about the self-force on a static electric charge was predicated 
on the horizon acting as conducting surface; therefore we should not be surprised that the 
force on the static scalar charge will be different. 



^ The factor l/u* is present even in a special relativity. Notice it is just the inverse of the Lorentz 
factor 7 = dt/dr = l/Vl — v"^. As a mnemonic, the strength of the "charge" depends on the 
spin of the field: for masses (spin 2) we have m = j^ni^Qst^ for electric charges (spin 1) we have 
e = 7*^ Crest) and clearly, by induction, for scalar charges (spin 0) we have q = 7~"^grest- 



in 



III. SELF-FORCE ON A STATIC SCALAR CHARGE 



A. local method 



We begin our calculation of the self- force with a picturesque description (Gedankenex- 
periment) of how such a measurement could be made. We imagine a test-charge with bare 
mass /iftare ^ud scalar charge q held fixed by a non-conducting system of mechanical struts 
outside the horizon of a Schwarzschild black hole. A non-conducting experimenter at the 
apex of a vertical, ballistic trajectory momentarily comes to rest with respect to the fixed 
charge. At this moment, she reaches out and measures the force required to hold the charge 
fixed, i.e. she measures the force needed to just lift the charge off the strut. The spacetime 
event B where/when the force is measured will be taken as the origin of the free-falling ob- 
server's coordinates x". (The over-bar denotes coordinates in the local, freely-falling frame 
of the observer [Q.) Clearly, by symmetry, we can choose our coordinate system such that 
the particle is located on the z-axis, and the freely-falling coordinate system is aligned so 
that the 2-axis coincides with the z-axis. Thus the only component of the force for the 
freely-falling observer to measure will be F^. Although this is an elaborate scheme to define 
the force measurement, working in the freely-falling frame where the charge is momentarily 
at rest is the surest way to establish unambiguously how the scalar field of the particle 
contributes to the renormalized mass. 

As no scalar charges have ever been observed, nor scalar fields measured, our assumption 
that the scalar field does not interact with the experimental apparatus (the struts) and the 
experimenter seems to be quite plausible. 

In order to compute the force in the free-falling frame of the observer, we need to be able 
to convert quantities from the isotropic coordinates of Eq. (|2.16|) to the coordinates of the 
observer. The defining feature of the free-falling frame is that spacetime is locally flat: i.e., 



In particular, the Christoffel symbols, r|- = {l/2)g°"^[g^;yp + g^p^^y — Qp^y^a] ~ 0[x'^] near 
the event B . Fortunately, SW have done the dirty work in finding the transformations from 
isotropic coordinates (un-barred) to the free-falling coordinates (barred). They are given 



(i/a^)^ = diag(-l, 1, 1, 1) = r/^^ 
9^-p{x^)=V^-p + 0[{x^f]. 



(3.1a) 
(3.1b) 
(3.1c) 



by 




(3.2b) 



(3.2a) 



These can be used to find 



1 



1 



1 



{l + M/2hf 1 



M 1 



z + 0[(x")2] , (3.2c) 
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where X(f) is the position of the charge as viewed in the freely-faUing frame 



(3.2d) 



and ag is the acceleration of the fixed charge as measured in the freely falling frame at the 
moment the experimenter comes to rest at the apex of her geodesic trajectory 



M 



1 



ag= 



62 (1 + M/2&)3(1 - M/2h) 



(3.2e) 
(3.2f) 



In the second line we have converted to Schwarzschild coordinates. Notice this is the same 



as the acceleration appearing in Eq. ( |1.2| ). We can use Eq. (|3.2|) to evaluate the scalar field 
and its derivatives in the free-falling frame 



\/(t,x) = - 
V{t = 0,x) 



|x-X(t)| 
■1 



1 - -agz + 0[{x^y] 



-q -n^ + 0\r X even number of (n') 'si 

[r 2 

Vt{t = 0, x) = — g 0[f° X odd number of (n^) 's] , 



V-,{i=0,5t) = -q 
V;iT(t = 0,x) = -g 



+ ^{n^n^ - S^^) + n^O[r° x even number of (rJ) 's] 
- — + 0[f^^ X even number of (n') 's] 



2 f 



r 



n^Olr" X odd number of (ra') 's 



(3.3a) 

(3.3b) 
(3.3c) 
(3.3d) 

(3.3e) 
(3.3f) 



where f = |x| and = x^/f. The spatial indices can be freely raised and lowered with 5*-'. 

Notice the extinction factor present in Eq. ( p. 16 ) does not extinguish the charge in the 
freely-falling frame Eq. ( |3.3| a), that is, in the freely-falling frame the dominant behavior of the 
field is simply (charge/ distance), independent of how deep the charge is in the Schwarzschild 
potential. 

We will compute the force required to hold the charge fixed by integrating the force 
density mm 



f = T 



'Zl3 



;/9 



(3.4) 



over the physical extent of the charged body at the instant of time {i = 0) when the 
measurement is made. More precisely, since we are describing the particle as a Dirac 6- 
function, we will integrate over an infinitesimal sphere of radius e centered on the particle 
and take the limit as e — >■ 0. The stress-energy tensor T^^ will have contributions from the 
bare mass of the particle and the scalar field; thus we have 



(strut) 



lim / [f^lf-o d X = lim 



1 


rpZ(3 


d^x 


Jf<e 







lim 



rpZl3 

(bare) ;/3 



d X + lim 

4=0 e^Ojf<f 



rpZf3 

(SF) ;^ 



d-^x 



t=0 



(3.5a) 
(3.5b) 



19 



The first term in the Eq. ( |3.5b| ) shows that the strut must support the bare weight of the 
particle. This term is present whether or not the particle is charged. Using the stress-energy 
tensor for the bare mass of a point particle located at b = 6z 



T, 



_ h'^{h)h^{h) ^3 
(bare) " Mbare ^ — 

g{b)u\h) 



(3.6) 



SW found the necessary force the strut must supply to support the bare weight of the particle 
to be 



(bare) 



lim 



T, 



zl3 



(bare) ;/3 



bare 



1 



1 



t=0 
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/^bare'^g 



(3.7) 



The second term in Eq. (|3.5b| ) involves the stress-energy tensor of the scalar field. In 
evaluating this integral we make use of Eq. ( |3.1| ) and note that the connection coefficients 
in this frame are 0[x"]; thus we can write 



F.opx= lim 

(SF) ,-^0 



lim< 

e^O lJf<i 



rpZl3 

(SF) ;;3 



t=0 



-(SF) 



t=o 



(3.8) 



We will denote the three contributions to Eq. ( |3.8| ) as -F(sfi)5 -^(sf2) -^(sf3) respectively. 
We also make use of Eq. ( p.l|) in writing the stress-energy tensor for the scalar field 



1 

An 



(3.9a) 
(3.9b) 



We now substitute Eq. (|3.9b|) into Eq. (|3^ ) and treat the terms in reverse order. The third 
term in Eq. (|3.8| ) gives a contribution of the form 



1 

47r 

-— lim 



arpzt ] 

(SF)J 



d-^x 



t=o 



o[x''v^v^] + o[ixyv^v, 



d-^x . 



(3.10) 



Using Eq. ( ^731) , we see the most singular terms come from (V^sV^f); therefore we have 



(SF3) 



lim 



— X [term with odd number of {n ) s] 

r 

-|- f° X [term with even number of (?t,')'s 



drdVL . 



(3.11) 



The first term contains an odd number of unit vectors, and therefore will vanish when we 
integrate over the solid angle. The second term vanishes as e — > 0, and similarly for higher 
powers of f. These two tricks are used repeatedly in evaluating the remaining integrals in 
Eq. (|3.8|). Thus we have 



'(SF3) 







(3.12) 



We now evaluate the second term in Eq. ( |3.8| ) using Eq. ( p.3| ) and ruthlessly discarding 
terms that do not survive the limit or the angular integration: 



-^(SF2)- 



-lim— / [VgV-Jf'dfdQ 

— q'^tta lim / —df q> n^n^dVt 
47r ^s^oJo r2 / 

-q^a„ lim / —df . 
3 ^e^oJo r2 



(3.13a) 
(3.13b) 
(3.13c) 



Using the divergence theorem and Eq. (|3.9b| ), the first term in Eq. (|3.8|) can be written 
1 



= — lim i 

An f-^o Jf=i 



rpzk 

^ (SF) ,fc 



n^f'^dn 



d^x = lim f 

t=0 e^Ojf=eL ~ 'Jf=0 



rpzk 

^(SF) 



f^n^d^l . (3.14) 



i=0 



All the terms except the first vanish either by / . . . dVt = 0, or lime_,o . . . = 0. Using Eq. ( p.3| 
and integrating the first term in Eq. ( p.l4| ) over the solid angle, we are left with 



(SFl) 



1 



Combining the results in Eqs. ( p.5b| ), ( |3.7] ), ( |3.12| ), ( p.l3cp and ( |3.15| ) we have 



(3.15) 



(strut) 



/"bare + — lira - + / —dr 



3 ^--0 Ve 

^2 







rco I 
/ibare + V / 

3 JO . 



9 ' 



(3.16a) 

(3.16b) 
(3.16c) 



where we have defined the leading factor in the bracket as the "renormalized" mass ^ren- The 
renormalized mass has the same functional form as the renormalized mass for the electric 



charge in the SW calculation. Converting to Schwarzschild coordinates [|14|, we get 



^ (strut) 



hi 



2M\-i/2 



+ { nothing depending on g } . 



(3.17) 



There is no self-force of the form seen in the second term of Eq. (|1.2| 



B. global method 

We verify our no-self-force result by means of a global, energy-conservation calculation. 
Suppose, instead of measuring the force on the charge while the charge is in place (as we did 
in the last subsection), the free-falling observer lowers the charge a small amount 5h. The 
work done on the experimenter will be 



14 



5W = -F'Sb = + M/2hf5h , 



(3.18) 



where we have used Eq. ( p.3|) to convert the free-falhng displacement 6b to an isotropic 
coordinate displacement 6b. The experimenter then converts this energy into a photon and 
fires the photon to asymptotic infinity. The energy received at infinity will be red-shifted 

(5E,eceived = ^f^oM . (3. 19) 

By conservation of energy this change in the system will be manifested by a change in 
asymptotic mass Ai of the system 

6M= -<5^,eceivcd = [1 - {M/2bf]F~'6b . (3.20) 

Thus we have 

1 5M 



l-(M/26)2 6b 



(3.21) 



We now use the total mass variation law of Carter |3^, which shows how the asymptotic 
mass will differ between two situations where the gravitational and matter status of the 
spacetime is slightly altered. In our case we compute the difference in asymptotic mass 
before and after we make the small displacement of the charge. The relationship is given by 

5M - ^6 A = ^5 / G° ^V^d?x + / G^^h^,^gd?x . (3.22) 
87r 87r J lovr J 

Here k is the surface gravity of the black hole, A is the area of black hole, G^'^ is the Einstein 
tensor, and h^i, is the difference in the metric between the two configurations. In Eq. ( p.22|) 



we have neglected terms involving the spin of the black hole. Throughout this paper we have 
assumed that the metric is unperturbed by the presence of the charge; therefore the last 
term in Eq. ( p.22| ) vanishes. For a Schwarzschild black hole the area term in Eq. ( p.22| ) is 



just the change in the mass of the black hole. During our slow displacement we will assume 
that no matter or radiation goes down the hole, and therefore this term will vanish. (We 
revisit this point at the end of the section.) Using Einstein's equation to write G^q = Sn T°q 
all that remains of Eq. (|3.22| ) is 



6M= 6 J T\ 

= S I r° ,,e)o V^g d'x + 6j TO p)o d^x . (3.23) 
The first integral can be computed by Eq. (|3.6|). SW give the result 

/ T(Le)0 d'x = /iba.c( j~j|^^^^ ) • (3.24) 



^^bare — 



Noting that the metric is diagonal, the scalar field is static, and employing the definition of 
the stress tensor for a scalar field Eq. ( p.9a| ), we can write the second integral in Eq. (|3.23 ) 

as 



SF = 



-'(SF)O 



-g (fx = — 



in 



1 

8^ 



-g d?x 



(3.25) 



The first term in Eq. (|3.25|) can be converted to two surface integrals: one over tlie liorizon, 
tlie otlier over a spliere at r — oo. Tlie field and its derivatives are well behaved on the 
horizon, but g-^^^—g vanishes there; therefore the integral on the horizon vanishes. The other 
surface integral vanishes in the limit r oo. Using the original field equation Eq. ( |2.1| ), we 
have 



.j),k 



Any 



-9P 



Thus we can write 



h)d^X 



(3.26) 



(3.27) 



Unfortunately, the scalar field is divergent at the source point b; therefore we must renor- 
malize. We do so by modeling our source as a charged, spherical shell of radius e, i.e. 



53(x-b) 



-9 



lim 



6\\^-h\-e) 



(3.28) 



Performing the integration and using the metric to convert the radius of the ball to free- 
falling coordinates e = (1 + M/26)^e, we have 



SF 



1 - M/2b\ 



1 + 



r-r lim 

M/2bJ f-o 



lim — . 



(3.29) 



As expected the functional form is the same here as in Eq. (|3.24| ). 

Combining Eq. ( |3.24| ) and Eq. ( |3.29| ) defining the renormalized mass 



and using Eq. 



and Eq. 



= /ibare + - lim — 

2 e-^o e 

, we have 



(3.30) 



(3.31) 



This agrees exactly with our previous calculation Eq. (|3.16c|) : no finite part of the self-force. 

We close this section with a pedagogical comment on the global energy conservation 
method for computing the force. We note that the calculation was predicated on the as- 
sumption that the area term in Eq. ( |3.22| ) vanished, i.e., 



K 



— 6 A = 5(mass of the hole) = . 

Stt 



(3.32) 



In terms of modern black hole theory, this is a valid assumption: No particles were dropped 
into the hole. The charge was displaced slowly, so no transverse fields {i.e. radiation) 
heated the horizon. Therefore, the area remains unchanged. However, in Sec. [Ill A| we 



gave a primitive derivation of the force which did not appear to exphcitly rely on any 
sophisticated properties of black holes. (Implicitly, we did assume that the mass of the 
hole remained constant when the observer wiggled the charge to measure the force.) An 
interesting interpretation of the two force calculations is to accept the primitive derivation 
in Sec. [ill A| as the correct force. Then, when we evaluate the right hand side of Eq. ( ^.23[ ) 
and show that it is the same as our local (primitive) force calculation, we have verified that 
the area of the black hole did not change when we lowered the charge. This tells us, in spite 
of the fact that energy of the scalar-static field (or the electric field for an electric charge) 
extends clear down to the horizon, that none of the energy near the horizon is pushed across 
the horizon when we lower the charge. 



IV. CONCLUSIONS AND DISCUSSION 



We have developed a formalism for constructing the Hadamard elementary solution for 
the Schwarzschild-Helmholtz equation ( p.4|) . This formalism was chosen because it expresses 
the singular nature of field in the near proximity of the charge. In the case of a static charge, 
we are able to find a closed-form expression for the field. We have used this expression to 
show (after mass renormalization) there is no self-force on a static scalar charge outside a 
Schwarzschild black hole 

Although we have patterned our discussion after Copson ^6] and SW JTSf, there are 
alternative methods for computing the scalar field and the self-force. For example, in the 
static limit, Linet has used generalized axially symmetric potential theory (GASP) ||3^ 
to derive the Green's function for the scalar field. However, this type of construction has 
not been extended to the Helmholtz-type equation depicted in Eq. (|2.4| ). Lohiya has 
demonstrated a concise method for determining the force on a static electric charge. Lohiya's 
method also uses the Copson-Linet closed-form expression for the electrostatic potential; 
however it is unclear how to extend Lohiya's method to moving particles. The formalism 
we have developed can be extended to moving charges. The recursion relation, Eq. ( p. 131 ), 
can be integrated with 7^ 0. Although it may be hard to find a simple summation of the 
results as in Eq. ( 2.15bD , it is possible to obtain the Green's function to the first several 
orders in F. 

Now that we have computed the (absence of) forces acting on a static scalar charge, let 
us assemble what is known about all the forces on a static charge outside a Schwarzschild 
black hole. In order to express this, let us slightly change the thought experiment. We will 
give the test-charge a mass /x, an electric charge e, and a scalar charge q. We will support 
the charge on a strut as before, but, instead of measuring the force supplied by the strut at 
some moment t = 0, we kick the strut out from under the charge and find the instantaneous 
acceleration of the falling charge in harmonic coordinates. After the particle begins to 
move there will also be a radiation-reaction force, so we must make the measurement at the 
moment we remove the strut. The metric can be used to convert quantities from free-falling 
(proper) coordinates to harmonic coordinates W^- The result is 







r dti \ 





rh + m 



(strut) 



(4.1) 
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Using this to convert Eqs. ( |1.2| ) and ( |3.17| ), and expanding in the post-Newtonian quantity 
M/rhi we have 



dtl 



M 



-1+4 — -9 — +16 — + known terms 
\rh/ \rh/ \rhJ 



e2 r 

2 r 



^V6f-V + -('-V 



H 1 — I + known terms 

WhJ 2 \rhJ 



+ 



q_ 

M 

M 



all terms are known to be zero 



, M\ 87/M\2 

2 1 — — + unknown terms 

rhJ 4:\rhJ 



(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 
(4.2e) 



Here Vh denotes the radial position of the particle in harmonic coordinates. In line (|4.2| a) we 
have recovered the velocity independent terms of the geodesic equation of motion expressed 
in harmonic coordinates. Since this is just a Taylor expansion of the first term in Eq. ( p.. 21) , 
we know these terms to all orders in M/r^. Thankfully, the first three terms are in agreement 
with the second post-Newtonian equations of motions. (See e.g. ||39|j40[| .) Line (^]2|b) is just 
the expansion of the second term in Eq. ( |1.2D , and thus we know these terms to all orders 
in M/vh- Line ( [4.2| c) is the scalar-charge part, which we have shown to vanish for all orders 



in M/vh. For moving charges, there will very likely be non-zero terms. In line ( |4.2| d), only 
the first two terms are known from second post-Newtonian calculations. (See e.g. [^ ,^ . ) 
The question remains, can the unknown terms in line ([4.2| d) be obtained by methods similar 
to those used to find the electric and scalar forces, that is, by looking at the field (metric) 
perturbations produced by the mass of the test-particle? Obviously there are a number of 
conceptual issues to tackle in answering this question. For example, when the metric itself 
is the perturbed field, can we define a freely falling observer in the same way as we did in 
the scalar-charge force calculation? When solving for the metric perturbation, how do the 



stresses in the strut affect the solution? This is currently under vigorous investigation 
Line (|4.2| e) represents higher order effects, such as additional forces on the particle due to 
the change in the metric produced by the electric field of the particle. Such terms will also 
arise from a gauge change, say, r r + 0[/i]. These terms are clearly second order in the 
perturbations. 
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APPENDIX A: RECAP OF ELECTROSTATIC CHARGE IN SCHWARZSCHILD 

SPACETIME 



In this section we summarize the results of Copson and obtain the closed-form 
solution for the electrostatic potential using our Hadamard construction. Using isotropic 
coordinates and assuming the field is strictly static, Maxwell's equations for the electrostatic 
potential can be written 



C'^Ao,ij{x) + CMoj(x) = -167re 



where 



= h{r) 



d 
dr 



In 



(26 + M)3 
(1 + M/2rf- 



1 - M/2r 



X-' 



(Al) 



(A2) 



Equation ( [Al| ) is in the same form as Eq. (|2.4|) ; therefore we proceed using Eq. ( p.9|) and we 

get 



(ci) _ r /2r' + M\ 3/2 / 2r - M 

^0 --I^TTmJ I: 



.2r' - M 



(A3) 



The superscript (el) denotes that these are parts of the electrostatic solution. Using the 
recursion relation Eq. ( |2.13|) , we have 



uf\ry)-. 



r7-(cl)/ i\ "^rrCcl) 2 

>(r,r)= --U^ ^7 , 
U3 '{r,r)= —U^ '7 



Ut'\r,r'y- 



16 
45 

128 



(A4a) 
(A4b) 
(A4c) 
(A4d) 



where 7 is the same as in Eq. (|2.14| a). Once again the summation is elementary: 

?7f)(l + 27r) 



y/T VI + 17 



(A5) 



As in the scalar case, we integrate the elementary solution against the source and obtain a 
particular solution to electrostatic field Eq. (|Al|) 



^r(x,b) 



bhTh - cos e 



{bh + M){rh + M) y^r2 - 2rhbh cos O + bl^M^ sin^ 6 



(A6) 



where we have switched to harmonic coordinates This is Copson's 1928 solution 
"for the potential of an electron in the Schwarzschild field". As in the scalar case we must 
ask: does the particular solution satisfy the boundary conditions? It is well behaved at the 



1Q 



horizon, so there is no problem there. However, as r/^ ^ oo the potential does not give the 
correct value 

Ar\r, ^ oo) ~ -^-^ ^ - . (A7) 

The fact that the field does not behave as e/rh for large Vh suggests by Gauss's law that 
we have found a solution with some additional charge lying around. However, our solution 
satisfies the homogeneous (source-free) equation everywhere outside the horizon except at 
the source point where there is a charge e. We must conclude that we found a solution 
with some charge on the horizon. In order to fix the boundary condition, we need to add a 
monopolar solution of the homogeneous equation, that is, we need to add an image charge. 
It is easy to see what is needed, 

/■homog eM 

^° -(r, + M)(6, + M)' ^^^^ 



and check that this satisfies the homogeneous equation outside the horizon. Linet 
noticed the discrepancy in Eq. (|A7|) and added this piece to Copson's result. Combining 
the two pieces gives the final result: the electrostatic potential for a fixed charge outside a 
Schwarzschild black hole 

. / ,x e I bhTh — M'^ cos 9 ,A , . 

Ao(x, b) = <^ , + M} . (A9) 

{bh + M){rh + M) y y^r2 - 2rf,bh cos O + bl-M'' sin^ 9 J 

The potential A^^^^ in Eq. ( [A^ ) is the particular solution constructed directly from the 
Hadamard elementary solution. If the Hadamard potential A^^^^ is taken to be the actual 
potential and the force calculation is carried out (i.e., repeat the SW calculation using the 
method similar to Sec. |111 A|), the resulting force is zero. This means that although our 



construction of the Hadamard solution was strictly a local calculation, when we summed 
the series we found a solution with just enough charge on the horizon to cancel the repulsive 
force in Eq. (|1.2|) . This also means that the repulsive force that SW found for the electric 
charge is due solely to the part of the potential ^0°™°^ which is tacked on to satisfy the 
boundary conditions. In other words, the second term in Eq. ( |1.2| ) is simply the force 
produced by the image charge on the horizon, and the force can be computed from Eq. 
directly 



self = 



■rh=b. 



h 



e'M e'M ^^^^^ 



This gives a physical interpretation of the repulsive force; the charge outside the black hole 
is repelled by an image charge inside the horizon. 
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APPENDIX B: COMPARING CLOSED-FORM SOLUTIONS WITH SERIES 

SOLUTIONS 



Equating our closed-form solutions for the scalar-static and electrostatic fields with the 
conventional infinite series solutions, we can obtain some interesting summation formulas 
that do not appear in the standard references [P^ -|4^. 

For a fixed point source Eq. ( |2.1|) is easily solved by separation of variables. The angular 
dependence is expressed by Legendre polynomials, and the resulting radial equation is also 
Legendre's equation. Equating the series solution to our close-form solution Eq. ( |2.16| c) we 
have 



^(xh, bh) = -gW f" r ^ (Bla) 

\ bh + M L2^_ 2r^hh cos O + hl-M"^ sin^ 9 



_ q b,-M ^ Pi{r,/M)QMM) , if r^ < b^ ] . . 

Here the Pi and Qi are the Legendre functions. This summation is a special case of Equation 
28 in MacRobert 131. 



If the field point is located on the horizon (M = Vh < h) or on the axis (6 = 0) we can 
verify this formula using the standard summation formula [^^44| 



OO 1 

Y,{2n+l)Qn{x)Pn{y) = |x| > 1 and |x| > |y| . (B2) 

Applying this summation formula to Eq. (HI), we get 



Ffhorizon) = 5 ^ . (B3) 

^ ' bh-Mcos^Vbh + M ^ ' 

Clearly the horizon is not a surface of constant "potential". This is in contrast with the 
electrostatic case where the horizon is a surface of constant potential. [See Eq. (|A9|) .] On 
the axis of symmetry (i.e., 6* = 0, so Pi{cos9) = 1 for all values of /) we can also use Eq. ( [B2[ ) 
to sum Eq. (0) 



Mr ■ \ ^ /bh -M 

y(axis) = --^ ; — T\- — . (B4) 

^ ' |rh-bh|Vbh + M ^ ' 

In the electrostatic case, the series solution can be similarly obtained by separation of 
variables. The radial functions are derivatives of Legendre functions. Equating the analytic 
expression with the series solution gives an interesting summation formula 

e f bhVh - cos 6* 1 , , 

^ox,b =- — — ^ , +M} B5a 

{b, + M) (r, + M)\ _ 2r,b, cos 6 + bl ~ sin^ 6 ^ 



^ V 2/ + 1 [ (r, - M){b, - M)Pl{r,/M) Q[{bn/M) , if t^<K 

h Kl + 1) ^ \ (r, - M){bh - M)PI{bh/M) Q'lirn/M) , if rj, > bh 



(B5b) 



91 



Here it is understood when / = we make the replacement 



^^^^-j^P;{rh/M) M (when / = 0) . (B6) 

Notice that the horizon {vh = M) is a surface of constant potential. The closed-form 
expression [modulo the homogeneous piece in Eq. was computed in 1928 by Copson 

26| , who expanded the result in terms of radial functions and discovered the summation 



formula. The series result was rederived by Cohen and Wald p2|, and Hanni and Ruffini 



23| . The asymptotic form of the solution can be seen immediately from the closed-form 



result, or from the series by noting that Qq{x — > oo) '^1/x. We see that 



e 



Ao(r ^ oo) ^ - , (B7) 



r 



which is the correct behavior. Equation ( Pl| ) can also be obtained from Eq. ( |B5| ) by differ- 
entiating and using Legendre's equation. 
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